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ABSTRACT 

With a luminosity > 10^ L© and a mass-loss rate of ~ 2.10"'' Moyr"', the red supergiant VY CMa truly is a spectacular object. Because of its 
extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, 
provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have 
obtained spectroscopy of VY CMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 - 672 micron. 
The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 
13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low 
as ~1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope. 
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1. Introduction 

The oxygen-rich red supergiant VY CMa is one of the most ex- 
treme objects of its kind. Its large mass-loss rate results in a 
very dense circumstellar envelope and an enormous infrared ex- 
cess dHarwit et al.l 1200 lb . Infrared observations are invaluable 
for in-depth studies of the complex chemical and dynamical 
processes at work in its circumstellar envelope. By this means, 
VYCMa has b een shown to harbour a rich, water-dominate d 
chemistry (e.g. iNeufeldet aH 119991: iPolehampton et all l2010l) . 
The unprecedented wavelength coverage and spectral resolution 
of the Herschel instruments will undoubtedly push our under- 
standing of the mass-loss related phenomena a step further 

In this context, VYCMa has been observed as part of 
the Herschel Guaranteed Time Key Programme "Mass-loss of 
Evolved StarS", MESS (Groenewegen et al., in prep.). The 
spectra discussed hereafter have been ob tained with the PACS 
and SPIRE instruments of Herschel (Pil bratt et al.ll2010h . The 
instrume nts and thei r in- orbit performance are described in 
Poglitschilin (l20Tol) and lGriffinet"an ( l20Toh . The SPIRE cal- 
ibratio n methods and accuracy are outlined by ISwinvard et al.l 
(120101) . 



2. Observations and data reduction 

The PACS observations of VY CMa consist of 5 spectral energy 
distribution (SED) scans, acquired on 3 and 13 November 2009 
(see Table [T]for details). All observations were performed with 
a non-standard version of the PACS-SED observing mode. The 



difference with the standard mode resides in the wavelength cov- 
erage, as our SEDs cover the complete PACS wavelength range 
in the four bands, in two observations. All observations were 
"chop-nodded", and contained one single nodding cycle, and one 
single up-down scan in wavelength. 

The data were reduced with the nominal pipeline, except that 
we combined the spectra obtained at the two nod positions after 
rebinning only. The flux calibration is based on a constant detec- 
tor response and the wavelength ranges affected by light leaks 
have been excluded from the present analysis. 

The original design of the PACS SED observing mode in- 
cluded Nyquist spectral sampling. However, in the early version 
that was used for these observations, a few wavelength ranges 
were still slightly undersampled. Consequently, rebinning the 
spectra at the nominal resolution of the spectrometer returns a 
number of empty bins. For the modeling, we worked at the nom- 
inal resolution, and used linear interpolation at the wavelengths 
of the empty bins. For the line identifications, we have rebinned 
the spectrum to a slightly degraded resolution (half Nyquistfl 

The SPIRE spectrum of VY CMa was obtained on 
Operational Day (OD) 123 and consisted of 30 repetitions. 
Each repetition consisted of 1 forward and 1 reverse scan of 
the FTS. The total on-source integration time was 3996 sec. 
The unapodized spectral resolution was 1.4 GHz (0.048 cm '), 
after apodizati on (us ing extended Norton-Beer function 1.5; 
Na ylor & Tahi^l2007l) this is 2.1 GHz (0.07 cm'). 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



' The standard PACS-SED was modified on OD305, and the full 
Nyquist sampling has been implemented since then. 
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Fig. 1. The complete PACS + SPIRE spectrum of VY CMa. The various spectral bands are presented in order of increasing wave- 
length. For the sake of clarity, PACS band B2A is not included, as it matches the wavelength range of PACS-B3A. The ISO-LWS 
spectrum is overplotted for comparison (Polehampton et al. 2010). For PACS, the spectrum of the central spaxel is displayed, after 
application of the point source correction (version 2.0). The strongest spectral features of the main molecules are indicated. 



3. Results 

3.1. Observational results 

Figure[T] shows the full PACS and SPIRE spectrum of VY CMa. 
More than 930 molecular emission lines have been detected, 
95% of which we could already identify. More than 400 of these 
lines originate from water vapour, but many different molecules 
contribute to the spectra: the identified species and their isotopo 
logues are '^CO 
H'^'O (p-Hi^ 

OH, 28siO, 2'^SiO, ^°SiO, ;.i>5, rn^iN, 3>^2 , ni^, \^u\- n^v^ is 
possibly present as we l l, as r eported by iPhillips et al.l (Il992h 
and Polehampton et al.l (l2010l) . but due to the high density of 
lines, detailed modeling is needed to give a definitive answer 
(Matsuura et al. 2010, in prep). Both low-excitation (e.g., CO 



'"^CO, C'^O, C"*0, ortho-H'^O (o-H'^'O), para- 
2 O (p-Hf O), o-H'^O, o-Hi'^O, p-H^'^O, CN, CS,"sO, NH3 



SiS, HCN. 



SO7 , H7S, rCiiL H3O+ is 



J=4-3, at 55 K) and high-excitation (e.g., o-H20(168,9-159_6) at 
5500 K) are found, allowing us to study the molecular envelope 
from the outermost layers down to the dust formation region 
close to the star In this way, the thermodynamical structure of 
the complex envelope and the chemical processes in the enve- 
lope can be determined. 

The flux calibrations of the two instruments match within 
15%, the different spectral ranges of a given instrument agree 
within 5%. The fluxes measured by PACS and IRAS at 60 and 
lOOyum are similar within 5 and 15%. The agreement between 
SPIRE and ISO-LWS is within a few percent at 195//m, and 
the PACS and ISO-LWS spectra coincide to 20 + 10% over all 
wavelengths. Note that the ISO -LWS fluxes wer e found to be 
high in comparison with IRAS. iGrv et al] (l2003h reported that 
the LWS fluxes were 20% above those of IRAS at the flux level 
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Table 1. Journal of observations. The target coordinates are Q'2000: 07h22m58.3s, (J2000: -25°46'03.2" 
PACS observations were 2x2 rasters with a step of 4" along each of the instrument axes (Fig.|2]l. 



On Operational Day 183, the 
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Fig. 2. Line peak intensity, expressed in Jy, for the o-H20(7i 6- 
62,5) line at 66.09 /im, in all PACS observations. The dashed con- 
tours represent in white, the diffraction limited telescope beam at 
66yL(m (5" diameter), and in black, the instrumental PSF at half 
maximum. The PSF is significantly wider due to the spatial un- 
dersampling of the beam (the spaxels are 9.4"x9.4"). The outer 
white contour is at 50% of maximum H2O intensity. The emis- 
sion is not resolved. The shape visible at low flux corresponds 
to the PSF wings. The white circles indicate the centers of the 
spaxels during each of the SED scans. The black dot in the cen- 
ter marks the object coordinates. 



of VY CMa. On the other hand, the PACS flux in FiglUmight be 
slightly underestimated, due to the extension of the source. 

Indeed, although it appears unresolved at first sight, in lines 
as well as in continuum (Fig. |2]l, VYCMa can not be consid- 
ered as a point source for PACS. Its flux, integrated over the 
field, shows an excess of 15+3% in the blue and 10+5% in the 
red compared to the flux measured in the central spatial pixel 
(spaxel) after application of the point source correction, for the 
observations of OD 173 (shown in Fig. [Til. The same compari- 
son, carried out on the Neptune SED that was obtained on the 
same OD shows Neptune (~2" diameter during the observation) 
to be compatible with a point source, hence confirming the ex- 
tension of the emission around VY CMa. The ratio between cen- 
tral spaxel and integrated flux was computed for all raster po- 
sitions. All of these ratios show comparable or larger excesses, 
hence ruling out the possibility that, on OD 173, the satellite was 
pointing significantly out of the peak of emission. 

Additionally, these central to integrated flux ratios, as well 
as dedicated measurements in the spectra, reveal that the line to 
continuum ratios are constant to within 5% over the whole field 



of view, at all PACS wavelengths, the only exception being the 
[Cii] line at 157.74/im, which has a stronger integrated flux due 
to its interstellar origin. Finally, we created maps of the line ra- 
tios between high and low-excitation transitions of water, e.g. o- 
H2O (7 1 6-62 5) 66.09/im / 0-H2O (33 o-22 1) 66A3fim in the blue 
and 0-H2O (73,4-72,5) 166.82yL(m/o-H20(2i,2-lo,i) 179.53;umin 
the red part of the PACS SED. These maps show no structure 
either, supporting the interpretation that all the emission is lo- 
calized in the central part of our field. Hence, although detailed 
comparison with a high quality PSF is still to be performed, all 
PACS observations are consistent with an extension of 5 to 10", 
which is compatible with the central pa rt of the nebula as im aged 
in the visible, near and mid-infrared bv lSmith et al.l (1200 ih . 



40 3.2. Modeling the envelope: CO, H2O, HCN and SiO 



Due to the rich water spectrum, most of the emission lines 
are blended. In addition, many line transitions have high op- 
tical depth. Consequently, the fractional molecular abundance 
can only be retrieved from a proper solution of the radiative 
transfer equation coupled to the statistical equilibrium equations. 
Although we realize th at the env elope around VY CMa has a 
very complex structure JSmith et al...2009.) . we have assumed a 
spherically symmetric envelope. Using this assumption, the av- 
erage density over afl directions of '^CO, '^CO, HCN, o-H'''0, 
p-Hj^O, and ^'^SiO can be determined. A more detailed study 
will be presented in Matsuura et al. (in prep.). 

We have used the non-local thermodynamic equilibrium 
(non-LTE) r adiative transfer code GASTRoNOoM (Deci netafl 
2006, '2010), with line lists and collisional rates as specified by 
Decin et al.. (.2010.) . The kinetic temperature is calculated from 
the balance of cooling and heating, with the main contributions 
coming from adiabatic expansion, grain-gas collisions, and CO 
and H2O ro-vibrational transitions. The conservation of momen- 
tum delivers the velocity structure. The model parameters are 
summarized in Table |2l a comparison between observed data 
and theoretical predictions is shown in Fig. [3] Specifically, we 
derive a mean mass-loss rate of 1.5 x lO^'^Moyr"', a fractional 
abundance of SiO of 3 x 10 ^ of HCN we find 3 x 10'' and of 
0-H2O being 3 x 10 "*. We derive a '^CO/'^CO isotopic ratio of 
60, and an ortho-to-para water ratio (OPR) of 1.27:1. Due to the 
high water optical depth, the latter remains rather uncertain. 



4. Discussion 



HCN: IZiurvs et all (l2009l) derived a HCN fractional abundance 
[HCN/H2] of 1.2 X 10"^ in the spherical outflow and of 7.5 x 
10"^ in the so-called red and blue flow for a mass-loss rate 
of 2 x lO^^^Moyr"' and a distance of 1500pc, being similar 
to our results. This abundance is orders of magnitude higher 
than the t hermodynarnic eq uilibrium value (TE) [HCN/Htot] of 
6 X 10 " (lDuarietal.il 19991) . A similar high abundance of HCN 
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Fig. 3. The continuum-subtracted PACS spectrum (black) between 173 and 188.5;um. The main contributing molecules and isotopes 
are identified. Features not yet identified are indicated with a 'U' . The first modeling results are shown, with one color per molecular 
species. 

Table 2. Model parameters for VY CMa. Numbers in italics are This indicates that the isotopologue ratio is (roughly) constant 
input parameters that have been kept fixed at the given value. in the spherical wind. 
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has also been derived fo r oxygen-rich (O-rich) asymptotic gi - 
ant branch (AGB) stars dZiurvs et al.ll2009l: iDecin et alj|2010h . 
Photochemical predictions for the outer envelope of O-rich AGB 
stars (Willacv & Millar 1997) are clearly too low by at least an 
order of magnitude. Pulsation driven shock models can account 
for this high concentration of HCN, yielding a HCN fractional 
abundance in the inner win d (between 2.2 and SR^,) rangin g 
from 2 X IQ-'' to 1 x IQ-^ (Duar i et al.|[T999l:ICherchneffll2006l) . 

Si O: The derived SiO abundance is consistent with the re- 
sults of lZiurvs et alJ (I20p7|). Both inner wind shock-induced non- 
equilibrium chemistry (iDuari et al.lll999^ Cherchnefi"2006') and 
outer wind photochemical processes (jWillacy & Millar 1997) 
predict almost the same abundance value. 

H2O: IZubko et al.l (l2004l) derived a water abundance 
[H2O/H2] of 4 X 10-"* from fitting ISO SWS and LWS ob- 
servations, with an OPR of 1:1, similar to our results. Low 
OPR ratios are also derived for the hypergiant NML Cyg by 
[Zubkoetal. (2004 andfor the O-ric h AGB star WHya ( 1:1.3 
bv iZubko & EhtziS (l2000h and 1:1 bv lBarlow et"ai1(ll996l) ). An 
OPR close to 1 : 1 refl ects a cold spin temperature, around 15 K 
( lMummaetal.lll987h . It is thought that the OPR value should 
correspond to the physical conditions where water was formed, 
i.e. the warm and dense regions of the envelope where the chem- 
istry is under thermodynamic al equilibrium and controlled by 
three body reactions. Hence, we could expect an OPR value of 
3. The large opacity of the water lines very sensibly reduces 
our se nsitivity to the OPR value, as also noted by .Barlow et alj 
(Il996l) . so that, even using Hj^O, we can at present not rule out 
the aforementioned expected value. 

'^CO/'^CO: From low-excitation CO lines, iMuller et al" 



( 2007 ) derived also a '^CO/'^CO ratio of 60, while lMilamet al. 



(l2009l) obtained a '^C/'^C value of 46 for the spherical wind. 



5. Conclusion 

Through the detection of more than 900 emission lines of vari- 
ous chemical species in VY CMa, PACS and SPIRE have shown 
their diagnostic strength to unravel the chemically complex 
mass-loss processes at work in evolved stars. Specifically, we 
have derived the HCN, SiO and water fractional abundances of 
the circumstellar material around V Y CMa, which also shows a 
suprisingly low ortho-to-para water ratio, close to 1 . The origin 
of the latter is not yet understood, but the Herschel observations 
of many evolved stars in the MESS programme might shed light 
on this topic. 
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